Dynamic Scaling at the Zero-field 2D Superconducting Transition 
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Zero-field current-voltage [I-V) characteristics of a thin ("two-dimensional") Bi2Sr2CaCu208+a 
crystal are reported and analyzed in two ways. The "conventional" approach yields ambiguous 
results while a dynamical scaling analysis offers new insights into the Kosterlitz-Thouless-Berezinskii 
transition. The scaling theory predicts that the universal jump of the I-V exponent a should be 
between z -\- 1 and 1. A value oi z ~ 5.6 is obtained for the dynamical critical exponent, and is 
corroborated by data from other 2D superconductors. A simple dynamical model is presented to 
account for the results. 



Perhaps no transition is better known than the Kos- 
terlitz-Thoulcss-Berezinskii (KTB) transition, which 
occurs in systems in the universahty class of the two- 
dimensional (2D) XY model and which is marked by 
the unbinding of topological excitations known as vor- 
tex pairs. Since the original theory was derived in the 
early 1970's, the KTB transition has been mentioned in 
systems as diverse as superfluids, [|j superconductors, 
2D lattices, [Q ferromagnets, |^ and liquid crystals. 

Evidence for Kosterlitz-Thouless-Berezinskii (KTB) 
critical behavior in two-dimensional superconductors in 
zero-field was first reported l^in 1981 and has since been 
widely studied in thin films, 2D Josephson Junction 
arrays, pO[ and more recently in high-temperature super- 
conductor (HTSC) uhrathin films. |l|] Evidence for KTB 
behavior (but not necessarily a transition) also abounds 
for multilayers and HTSC single crystals. 

Two main approaches can be used to investigate elec- 
tronic transport phenomena in superconductors near the 
KTB transition. In the first, more "conventional" ap- 
proach, there are two principal signatures of KTB be- 
havior. |^,|l^ First, the ohmic resistance should have 
a unique temperature dependence above the super- 
conducting transition temperature, Tktb- R(T) oc 
exp[—2^b/{T/TKTB — 1)], where 6 is a non-universal 
constant. Second, current-voltage (I-V) isotherms 
should be described by a power law, V (x I", at low cur- 
rents. § Below Tj^TB, 01 decreases linearly with increas- 
ing temperature, makes a "universal jump" from three to 
one at T = Tj^tb, 11] then remains ohmic (a = 1) for 
all T > Tktb- 

There are several difficulties with this approach. For 
example, the jump in a is strictly defined only in the 
7^0 limit, making it difficult to detect experimentally. 
Indeed, many of the original papers [fzj-p^ do not ob- 
serve or report universal jumps. Furthermore, it cannot 
address the eventual upturn of ohmic isotherms, observed 
near Tktb- |1 (See Fig. 1 below, and Refs. [|||Jl|,|l|l.) 
Similar behavior is also observed in finite field near 



the superconducting transition in bulk superconductors, 
where it is attributed to the competition between ther- 
mal and current-induced effects. |jl^;Q In the second 
approach, in which one applies the techniques developed 
for such bulk systems, it is possible to resolve both of 
these issues. 

The second approach treats both low and high current 
behaviors, via dynamic scaling. The main ideas were 
first discussed for 3D XY and glassy critical phenomena 
in bulk superconductors, [l^ In this analysis, dynamic 
phenomena near the critical point involves a correlation 
time scale t cx which like the correlation length ^, 
diverges at Tktb- [z is the dynamic critical exponent.) 
The scaling approach has two distinct advantages over 
the conventional theory. First, it extends the transport 
theory to finite currents, and is able to describe the ob- 
served upturn of ohmic isotherms. Second, it admits dy- 
namics which are non-diffusive. This issue has recently 
received much attention in bulk HTSC's near the 3D XY 
critical point, where unexpected (non-diffusive) dynam- 
ics have been encountered. 

For 2D superconductors, the I-V curves should scale as 
V = I^-'x±iI^/T), ll|,0 where is the seal- 

ing function for temperatures above (below) Tktb- The 
asymptotic behaviors of x±{^) can be deduced from the 
conventional KTB theory: lunx^o x+ix) = const, (ohmic 
limit); lim^j^o X- (2^) exp [— aln(l/x)] (thermally acti- 
vated vortex pair unbinding j2l)); lim^j^oo X±(3^) oc 
(critical isotherm). The univeral jump is expressed as 
the difference in (log-log) slopes between the two asymp- 
totic limits of x+ (x) - The distinctions between the scal- 
ing and the conventional approaches are apparent: (i) 
R{T) oc exp[-Zy/b/{T/TKTB - 1)] (2 replaced by z); 
and (ii) the "universal" jump of a is now from z-\-l to 1, 
not from 3 to 1. Of course the two approaches are equiv- 
alent for diffusive dynamics (z — 2). Despite the recent 
widespread application of scaling techniques to TV char- 
acteristics near a critical transition, jist we are familiar 
with only two cases where these ideas have been applied 
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to 2D superconductors. |20 2^ In the former, the dy- 
namic universaUty class was not exphcitly studied while 
in the latter, the issue was treated only peripherally. 

In this Letter, we apply both approaches to I-V curves 
from a thin crystal of Bi2Sr2CaCu208+5 (BSCCO). It is 
found that the conventional approach gives rather am- 
biguous results. However, by applying dynamical scal- 
ing, with z as a fitting parameter, the data can be made 
to collapse, obtaining z ~ 5.6 ± 0.3 {not z = 2). To 
check the universality of our results, we also apply the 
scaling analysis to other I-V data sets obtained from the 
literature, including both HTSC and conventional 2D su- 
perconductors, with an average result of z ~ 5.7. 

BSCCO crystals were prepared using a standard self 
flux method. A single crystal was selected and cleaved 
using a Scotch tape technique. ||2^ Thin 1500 A crystals 
were obtained with atomically smooth surface areas of at 
least (50 micron)^. Q Thicknesses were measured using 
an atomic force microscope. Four silver pads were pho- 
tolithographically patterned on top of the crystal. The 
crystals were postannealed at 600 C in 1.5 scfh flowing 
O2 for 45 minutes to achieve optimum oxygenation. Gold 
wires were ultrasonically bonded to the silver pads where 
they extended beyond the crystal. 

R{T) and I-V characteristics were measured using a 
standard four probe inline geometry and a low frequency 
lock-in technique (16.9 Hz). A temperature controller 
was used to insure a temperature stability of better than 
10 mK. Although current was injected only through the 
top crystal surface, the thickness of the sample assured 
uniform current distribution. Data was taken on three 
different samples. Data presented here is from a crystal 
of dimensions 0.2mm x 0.5mm x (1500 ± 500 A). For 
thicker samples, KTB behavior could not be observed. 

Fig. H shows zero field {H ^ 100 mGauss) I-V 
isotherms, corresponding to 100 mK steps, from T — 80.1 
K (left side) to 76.9 K (right side). The curves dis- 
play features typical of 2D samples. The high temper- 
ature (> Tktb) isotherms exhibit positive curvature 
as they cross over from ohmic to non-ohmic behavior, 
while the low temperature (< Tktb) isotherms exhibit 
negative curvature. The critical isotherm separates the 
two regions at T = Tktb, and describes a strict power 
law: V oc 7^+^. In Fig, [l, and for other published 2D 
isotherms {e.g., Refs. [plJTq^ 1^), identification of the crit- 
ical isotherm shows that z is clearly larger than 2. 

We first perform a conventional analysis. The ohmic 
resistance data R{T) is fit to the KTB prediction, with 
apparent success, obtaining Tktb — 78.87 K and b — 1.3. 
The TV exponent a{T) is then determined by fitting 
F oc /" to the lower portion (2 to 10 nV) of each 
isotherm. In Fig. ^ high temperature (low current) 
ohmic behavior, expected to persist for all T > Tktb, 
is obscured by noise already for T < 79.7 K. Such ex- 
perimental limitations lead to significant errors in a{T) 
near and below Tktb, thus smearing out the universal 



jump until it becomes unrecognizable, as in Fig. |l] (in- 
set). These same difficulties affect R{T) measurements, 
reducing the accuracy of the obtained fitting parameters. 




10-6 L 

100 1000 
I{HA) 

FIG. 1. TV isotherms from a thin Bi2Sr2CaCu208+5 crys- 
tal for temperatures 76.9 K to 80.1 K. [INSET: a{T) deter- 
mined by a fit to oc I" in the voltage range 2 to 10 nV.] 

In the dynamic scaling analysis, we are not restricted 
to low currents, so that voltage sensitivity no longer be- 
comes an issue. It is convenient to rewrite the scaling 
ansatz as [M 
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where £±(2:) = x/x±{x^^^)- Note that above Tktb, the 
scaling variable can also be written as a; = I^ / R{T)T^ . 
To proceed with the analysis, C(T) must be speci- 
fied. We will assume that the KTB form, ^(T) oc 
exp[^yb/ {T/Tktb — 1)], provides the most efficient pa- 
rameterization of the correlation length. We further as- 
sume that ^(T) is symmetric about T — Tktb (modulo 
some prefactor). The fitting parameters for Eq. (|l|) 
then become z (universal), and Tktb and b (nonuniver- 
sal). The following three requirements must be fulfilled 
self-consistently in our scaling procedure: (i) V oc I^^^, 
along the critical isotherm T — Tktb', (h) R{T) oc 
in the high temperature range (^ 79.7 K), where ohmic 
R{T) can be obtained; and (iii) scaling collapse of the 
TV isotherms, according to Eq. (|l|). The first two points 
place tight constraints on the fitting parameters. Within 
these constraints, the parameters are tuned to satisfy the 
final point. 

The outcome of such a fitting procedure for the data 
of Fig. |l| gives the following results: Tktb = 78.87K, 
b ~ 0.57, and z = 5.6 ± 0.3. The scaling collapse, shown 
in Fig. (H), is convincing in both the upper (T < Tktb) 
and lower (T > Tktb) branches. For many isotherms in 
Fig. § the noisy (low current) portions of the curves do 
not scale; this may be due to electromagnetic screening 
(finite size effects). |15 However the same difficulty is not 
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incurred for the other data analyzed below. (See Fig. |^.) 
The caution expressed earlier for the low temperature 
branch ||2l| does not seem to affect the scaling collapse. 
The inset of Fig. || shows the results of a similar two- 
parameter fit to the data, obtained by setting z — 2. 
The results are extremely poor, indicating that z> 2. 
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FIG. 2. The I-V curves of Fig. § scaled with Eq. [IN- 
SET: The same scaling for z — 2.] 

To check our results, we have applied the same scal- 
ing procedure to other data from the literature. The 
first data set ||l^ is from a YBCO/PrBa2Cu307-5 mul- 
tilayer in which the YBCO layers have a thickness of 
24 A and are electrically isolated from one another by 
PrBa2Cu307_5 barrier layers. The scaling procedure 
now leads to the results Tktb = 32.0 K, b = 39.3, and 
z = 5.6 ± 0.3, with the I-V collapse shown in Fig. ||a. 
Analysis of I-V data Q from a 12 1 YBCO monolayer 
gives Tktb = 18.5 K, & = 22.29, and z = 5.8 ± 0.3 and 
also collapses the data. (See Fig. |^b.) The final data set 
corresponds to a conventional 2D superconducting sam- 
ple j|] (Hg-Xe alloy) and collapses for Tktb — 3.04 K, 
b = 9.64, and z = 5.6 ± 0.3 as shown in Fig. |[:. The 
apparent agreement between these four diverse samples 
suggests that the mean result, z ~ 5.7, may indeed be 
universal for superconductors. However, it was not pos- 
sible to probe the universality of the scaling functions, 
because of the non-universal parameter b appearing in ^. 

The allure of the KTB transition lies in the simplic- 
ity and integrity of the pair-dissociation paradigm. It 
is therefore unsettling that the corresponding dynamic 
paradigm could fail, as signaled above through the obser- 
vation of z <^ 2. To better understand this situation, it is 
helpful to reconcile the intuitive ( "conventional" ) picture 
with the scaling description. The conventional paradigm 
intrinsically involves the presence of dissociated or "free" 
vortices, with density Uf. In a small but finite current, 
free vortices are independently driven across the sample, 
inducing dissipation consistent with the Bardeen-Stephen 



theory: R oc Uf. In larger currents, non-ohmic effects 
may arise from steady state dissociation-recombination 
processes, although dissipation is still understood in 
terms of driven, non-interacting vortices. To make con- 
nection with the scaling approach, we compare the di- 
mensional statement nf cc C '^, which cannot be affected 
by dynamical processes, with the 2D ansatz Q R oc 

— z/2 

This leads immediately to the relation R oc , which 
reproduces the conventional relation when z — 2. Us- 
ing this relationship in the derivation of the conventional 
signatures, one finds R{T) oc exp[— z y^b/ {T /Tktb ~ 1)] 
and that a jumps from z -|- 1 to 1 at Tktb, consistent 
with the dynamic scaling results. 
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FIG. 3. The I-V curves of (a) Ref. [@ (b) Ref. || and 
(c) Ref. Q scaled with Eq. (^). Data sets (b) and (c) have 
been shifted arbitrarily. 

The description in terms of a general z now allows for 
other more sophisticated dissipation phenomena. In par- 
ticular, the observed value of z ~ 5.7 cannot be mistaken 
to be diffusive, and suggests transport processes which 
are highly collaborative in nature. We speculate that 
for superconductors, driven vortices may never be "free" 
in the conventional sense. Rather, the dominant trans- 
port mechanism may involve collaborative dissociation, 
in which pairs of vortices can unbind only by exchanging 
members with neighboring pairs. In contrast with the 
conventional picture, which involves only the two mem- 
bers of the dissociating vortex pair, the collaborative dis- 
sociation scenario involves four "mobile" (as opposed to 
"free") vortices: the dissociating pair, and the two neigh- 
boring vortices with which they recombine. Assuming a 
steady state recombination process analogous to the con- 
ventional one, 1^ we now obtain rif ~ nf/r oc n'j, or 
r oc with a dynamic exponent of 6 rather than 2. We 
emphasize that the relatively large observed value of z 
should not be attributed to pinning effects, as the phe- 
nomenon is observed in a variety of superconductors and 



3 



must be regarded as intrinsic. 

Finally, we contrast the conclusions of this work with 
Ref. ||l5|, where it was suggested that the upturn of 
isotherms near Tktb, from ohmic to non-ohmic be- 
havior, may be attributed to electromagnetic screening, 
which produces free vortices even below Tktb- Such 
effects, which result from the 2D penetration depth be- 
coming smaller than the sample dimensions, should occur 
in some samples of finite thickness, ultimately destroy- 
ing the KTB transition. However, in analogy with bulk 
superconductors, the disputed behavior may instead re- 
flect only the non-ohmic behavior associated with finite 
currents. In this scenario, the upturning ohmic curves 
would correspond to T > Tktb, and the KTB transi- 
tion would remain intact. The non-ohmic behavior can 
be understood as follows. On either side of the tran- 
sition, external currents tend to dissociate large vortex 
pairs through (highly collaborative) activation processes. 
For T > TfCTB, many large pairs are already unbound, so 
that small currents, which tend to dissociate large pairs, 
have little effect. However, beyond some characteristic 
current scale (Jq ^ T/^), current-induced dissociation 
becomes the dominant dissipation mechanism, producing 
non-ohmic effects. It is worth noting that the disputed 
TV isotherms of Ref. are found to scale nicely with 
the techniques used here. (See Fig. ^3.) 

In summary, a conventional analysis of electronic 
transport data is contrasted with a dynamical scaling 
analysis for several types of 2D superconductors. The 
latter approach is designed to include non-ohmic effects 
associated with finite currents. The conventional analysis 
obtains inconclusive results while the dynamical scaling 
analysis yields a collapse of the I-V data. The results 
show that dynamical critical behavior in 2D supercon- 
ductors may be more interesting than previously antici- 
pated. The scaling generalization of the universal jump 
of the TV exponent a is from a value of z -I- 1 to 1. Our 
analysis gives a dynamical exponent of z ~ 5.7, rather 
than the expected value of z = 2, while a simple dynam- 
ical model predicts z ~ 6. Our results may pertain more 
directly to superconductors, where previous reports of 
anomalous vortex diffusion have been made, [^7|] rather 
than to pure XY models, where there is no indication of 
non-diffusive dynamics. It remains to be seen how the 
present results are refiected in 2D Josephson junction ar- 
rays and granular films. 
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